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ABSTRACT 

The redshifted ultraviolet light from early stars at z ~ 10 contributes to the 
cosmic near infrared background. We present detailed calculations of its spectrum 
with various assumptions about metallicity and mass spectrum of early stars. 
We show that if the near infrared background has a stellar origin, metal-free 
stars are not the only explanation of the excess near infrared background; stars 
with metals (e.g. Z = 1/50Z Q ) can produce the same amount of background 
intensity as the metal-free stars. We quantitatively show that the predicted 
average intensity at 1-2 /im is essentially determined by the efficiency of nuclear 
burning in stars, which is not very sensitive to metallicity. We predict vl v j ~ 
4 — 8 nW m -2 sr™ 1 , where p* is the mean star formation rate at z = 7 — 15 
(in units of M Q yr _1 Mpc -3 ) for stars more massive than 5 M Q . On the other 
hand, since we have very little knowledge about the form of mass spectrum of 
early stars, uncertainty in the average intensity due to the mass spectrum could 
be large. An accurate determination of the near infrared background allows us 
to probe formation history of early stars, which is difficult to constrain by other 
means. While the star formation rate at z = 7 — 15 inferred from the current 
data is significantly higher than the local rate at z < 5, it does not rule out the 
stellar origin of the cosmic near infrared background. In addition, we show that 
a reasonable initial mass function, coupled with this star formation rate, does 
not over-produce metals in the universe in most cases, and may produce as little 
as less than 1 % of the metals observed in the universe today. 

Subject headings: cosmology: theory — diffuse radiation — infrared: galaxies 
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1. INTRODUCTION 

When and how was the universe reionized? These questions have actively been studied 
almost purely by theoretical means (see Barkana & Loeb (2001); Bromm & Larson (2004); 
Ciardi & Ferrara (2005) for recent reviews), as currently there are only a very few observa- 
tional probes of the epoch of reionization: the Gunn-Peterson test (Gunn & Peterson 1965; 
Becker et al. 2001), polarized light of the cosmic microwave background on large angular 
scales (Zaldarriaga 1997; Kaplinghat et al. 2003), which has been detected by the Wilkinson 
Microwave Anisotropy Probe (WMAP) (Kogut et al. 2003), and temperature of the inter- 
galactic medium (Hui & Haiman 2003). Over the next decades, we are hoping to detect 
the first sources of light directly with the next generation of space telescopes, such as the 
James-Webb Space Telescope (JWST). More ambitious are mapping observations and mea- 
surements of the power spectrum of fluctuations of the 21-cm line background from neutral 
hydrogen atoms during reionization (Ciardi & Madau 2003; Furlanetto, Sokasian, & Hern- 
quist 2004) or even prior to reionization (Scott & Rees 1990; Madau, Meiksin, & Rees 1997; 
Tozzi et al. 2000; Iliev et al. 2002), which offer very powerful probes of detailed history of 
the cosmic reionization. 

It has been pointed out that the mean intensity (Santos, Bromm, & Kamionkowski 
2002; Salvaterra & Ferrara 2003; Cooray & Yoshida 2004; Madau & Silk 2005) as well as 
fluctuations (Magliocchetti, Salvaterra, & Ferrara 2003; Kashlinsky et al. 2004; Cooray et 
al. 2004) of the near infrared background potentially offer yet another window to the epoch 
of reionization. The logic is very simple: suppose that most of reionization occurred at, say, 
z = 9. The ultraviolet photons (A ~ 1000 A) produced at z = 9 during reionization will then 
be redshifted to the near infrared regime (A ~ 1 jum). In other words, a fraction of the near 
infrared background (whether or not observable) must come from the epoch of reionization, 
and there is no question about the existence of the signal. (Of course the existence of 
the signal does not immediately imply that the signal is actually significant.) It is therefore 
extremely important to understand the near infrared background in the context of redshifted 
UV photons and examine to what extent it is relevant to and useful for understanding the 
physics of cosmic reionization. 

Has it been detected? All of the theoretical proposals were essentially motivated by 
the current measurements of the near infrared background, which suggests the existence of 
an isotropic background after subtraction of the zodiacal emission (Dwek & Arendt 1998; 
Gorjian, Wright k Chary 2000; Wright & Reese 2000; Wright 2001; Cambresy et al. 2001; 
Matsumoto et al. 2005; Kashlinksy 2005). Since the zodiacal emission is ~ 3 times as large 
as the inferred isotropic component, one should generally be careful when interpreting the 
data. Although the inferred background is isotropic at the first order, significant fluctuations 
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still remain at the level of AI/I ~ 1/4 (Kashlinsky & Odenwald 2000; Kashlinsky et al. 2002; 
Matsumoto et al. 2005), which requires further explanations. The most intriguing feature of 
the current observational data is that the inferred background seems too large to be accounted 
for by the integrated light from galaxies (Totani et al. (2001), see also Figure 12 of Matsumoto 
et al. (2005) and references therein for the compilation of the galactic contribution). It is 
thus tempting to speculate that the bulk of the near infrared background (aside from the 
zodiacal light) actually comes from stellar sources at the epoch of reionization. 

In this paper, we carefully examine the near infrared background from early stars. While 
our approach is similar to that of Santos, Bromm, & Kamionkowski (2002), which has been 
adopted by most of the subsequent work, our goal is to (1) simplify physics and improve 
calculations, (2) explore different metallicity and initial mass spectra, (3) understand robust- 
ness of theoretical predictions, and (4) provide a simple relation between the cosmic near 
infrared background and star formation rate. The focus of this paper is the mean intensity: 
we will discuss fluctuations in the forthcoming paper (Fernandez et al., in preparation). This 
paper is organized as follows. In § 2, we develop the basic formalism and summarize rele- 
vant emission processes such as stellar emission and reprocessed light, the latter including 
Lyman-a, two-photon, free-free and free-bound emission. In § 3, we examine energy spectra 
of various emission processes from early stars. In particular we explore differences in the 
energy spectrum between various assumptions about metallicity and initial mass spectrum 
of early stars. In § 4, we calculate the spectrum of the cosmic near infrared background. In 
§ 5, we compare the prediction to the current observational data and discuss implications 
for the star formation rate at z = 7 — 15. In § 6, constraints from metal production from 
the first stars are calculated. In § 7 we discuss the other constraints from a collapse fraction 
of dark matter halos. We conclude in § 8. 

2. STELLAR EMISSION AND REPROCESSED LIGHT 
2.1. Basic Formalism 

We calculate the background intensity, I U1 as (Peacock 1999) 



where v is an observed frequency (which is in the near infrared band: for A = 3 — 1 /xm, 
say, v — 100 - 300 THz or his = 0.414 - 1.24 eV), H(z) is the expansion rate at redshift 
z [dtjdz = — [H(z)(l + z)}" 1 ), and p(u, z) is the volume emissivity in units of energy per 
unit time, unit frequency and unit comoving volume. There are several contributions to 
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the emissivity. One is the continuum emission from stars themselves, p*, which is nearly 
a black body spectrum, and the others are reprocessed light of ionizing radiation: a star 
ionizes neutral gas in its neighborhood and a series of recombination lines, p hne ) emerge. 
The ionized gas (or nebula) also emits free- free and free-bound continuum emission, p cont ) 
as well as two-photon emission, p 21 . In Appendix A, we derive the formula for the volume 
emissivity as (Eq. [A9]) 



and f(m) is a mass spectrum (specified later in § 3.1; for the precise definition, see Ap- 
pendix A), and is the mean stellar mass (Eq. [A5]). It is important to note that this 
formula has been derived assuming that the stellar main-sequence lifetime, r(m), is shorter 
than the Hubble time, and corrected for dead stars which do not contribute to the volume 
emissivity. Here, L,"(m) is a time-averaged luminosity (over the main-sequence lifetime) in 
frequency interval [u, v + dv\ for a radiative process of a, and (e") is the key dimensionless 
quantity which represents a ratio of the mass-weighted average 1 of total radiative energy 
(including stellar emission and reprocessed light) to the stellar rest mass energy, in unit fre- 
quency interval. In other words, (e°) represents the mass-weighted mean radiative efficiency 
of stars. This formulation is useful as one can immediately see that each contribution is 
simply given by the star formation rate (which depends on z) and a typical radiative effi- 
ciency (which does not depend on z). While p* is very uncertain and will be constrained by 
a comparison to the observational data, one can calculate (e°) robustly for a given popula- 
tion of stars using simple physical arguments. It will be shown in the subsequent sections 
that f(e") is always of order 1CT 3 when averaged over the main-sequence lifetime, which 
can be understood with simple energetics. Initially, energy must be generated by nuclear 
burning in stars. While the rest mass energy of 1 M is as big as 1.8 x 10 54 erg, only a 
fraction will go into radiative energy. For example, in the Sun only 0.07% of the rest mass 
energy is converted to radiative energy over its main-sequence lifetime. The nuclear burning 
efficiency 2 depends on stellar mass only weakly at large masses. Our detailed calculations 
below confirm this simple argument, and thus the uncertainty in the predicted amplitude of 
radiative efficiency is small for a given mass spectrum of stars. 



throughout this paper, we shall use (} to denote the mass-weighted average. 

2 By "nuclear burning efficiency", we mean the bolomctric energy of stellar emission before absorption 
per stellar rest mass energy, L bol T/mc 2 . 
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Using the expected radiative efficiency of stars, we obtain 

dy 10 2 [ My) 



vl v = 13.1 nW nT 2 sr^ 1 
where y = 1 + z and 



y 2 £(y) 



M Q yr- 1 Mpc- 



E 



io- 



£{y) = hy/n m y 3 + fi A ~ 11.83 ( 



0.14 



io; 



(4) 



(5) 



for the redshift range of interest. Thus, without any detailed calculations, one can show that 
the cosmic near infrared background from early stars at z ~ 10 should be approximately 
given by vl v j p* ~ 10 nW m~ 2 sr _1 where p\ is in units of M & yr -1 Mpc~ 3 . 



2.2. Stellar Contribution 

To simplify calculations, we assume that the stellar spectrum is a black body with the 
Lyman continuum photons completely absorbed: 

T(m\ - ! ^Rl(jn)B v [T eff {m)i hv < 13.6 eV 

LAm} ~ \ 0, hv> 13.6 eV (6) 

where is a stellar radius and T e ff is the effective temperature, and B u [T e ff(m)] is a black 
body spectrum given by 



2hu 3 / 



,2 



B„[T e// (m)] = MW) _ V (7) 

Note that the stellar spectrum (before absorption) above 13.6 eV (which determines the 
number of hydrogen-ionizing photons) is significantly different from a black body; thus, we 
do not use a black-body spectrum to calculate the number of ionizing photons, but use more 
detailed calculations by Schaerer (2002) (see § 3.2). The stellar spectrum just below 13.6 eV 
is also different from a black body because of a cluster of absorption lines of Lyman series; 
however, we ignore this effect and keep our calculations as simple as possible. (One can 
always use a more precise stellar spectrum for a better accuracy.) 



2.3. Free-free and Free-bound Contribution 



The free-free and free-bound continuum luminosity is given by 



rrV) = £vQH{m) , (8) 

n e n p a B 
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where Q H is a time-averaged production rate of hydrogen ionizing photons (the average 
number of ionizing photons produced per unit time), n e and n p are the number density of 
electrons and protons, respectively, as is the case-B recombination coefficient for hydrogen 
(Eq. [5-14] of Spitzer (1978) with n = 2 and Z = 1) given by 



a B 



2.06 x 10 

T l/2 



-11 



^(Tg) cm 3 s" 



(9) 



and <t>2{T g ) is a dimensionless function of temperature tabulated in Table 5.2 of Spitzer 
(1978). Here, T g denotes gas temperature in units of Kelvin. In principle, to calculate T g one 
has to equate the energy gain and loss to find out equilibrium temperature. While T g varies 
depending on stellar temperature (or hardness of a stellar spectrum which determines photo- 
heating), we shall assume T g = 20000 K regardless of the stellar temperature, which should be 
a good approximation for our purposes. For this temperature we find 2 (T g = 20000 K) ~ 1. 

The quantity Q H {rn)/{n e n p a.B) is the volume of the Stromgren sphere (see text below 
Eq. [12]), and e v is the total volume emissivity including free-free and free-bound emission 
(Eq. [6.22] of Dopita & Sutherland (2002)): 



e ~hv/kTg 

e v = ATin e n p ^ c erg cm" 3 s" 1 Hz" 1 , 

Tn 



(10) 



where 7 C is the continuum emission coefficient including free-free and free-bound emission: 



7c = fk 



9ff + J2 



n=2 



n 



-9fb(n) 



(11) 



where x n = Ry/ (kT g n 2 ), ~g ^ and gbf{n) are the Gaunt factors for free-free (which is thermally 
averaged) and free-bound emission, respectively, and /& is the collection of physical constants 
which has a numerical value of 5.44 x 10" 39 in cgs units. Note that we have ignored the 
helium contribution and assumed complete ionization for computing j c . As a free-bound 
transition to n — 1 will not be considered in the case-B recombination, the summation is 
taken from n — 2. (This is because all photons that recombine directly to n — 1 are strongly 
absorbed by neighboring hydrogen atoms and immediately ionize them.) We then obtain 



L c u ° n \m) ~ 3.32xl0 22 erg s" 1 Hz" 1 



10 49 s" 



Ry 

~ 9 " + W, 



™ e Ry/(kT g n^ 

2^ 9fb(n) 



9 n=2 



rr 



MT 9 )e- hu/kT °. 
(12) 



The continuum luminosity does not depend on the number density of electrons or protons. 
This is an immediate consequence of the Stromgren sphere: while the higher number density 
implies the larger emissivity, it also implies the larger recombination rate and the smaller 



-7- 



ionized region. These two effects cancel out exactly, making luminosity independent of the 
number density. Of course, this approximation breaks down in the intergalactic medium 
(outside of halos) in which ionization fronts do not fill the Stromgren sphere (Shapiro & 
Giroux 1987). Our calculation assuming the Stromgren sphere is accurate if the bulk of 
luminosity comes from nebulae around stars inside the host halos, while it should give a 
robust upper limit on free-free and free-bound luminosity otherwise. 

Finally, we need to compute the Gaunt factors. For the parameter space we are inter- 
ested in, 

- - hJ W-°^ ™ 
7 2 = §- g = °m> (14) 

both Gaunt factors are approximately constant and given by (Karzas & Latter 1961) 

9 ff - LI, (15) 
g fb (n) ~ 1.05, (16) 

which are accurate to within 10%. 



2.4. Line Contribution 

The line luminosity is given by 

L l T\ m ) = hv ^ v ~ Vi)ni(m), (17) 

i 

where 4>i{y ~~ v i) is the line profile and n« is a photon production rate at a line i. Since 
the intergalactic medium is optically thick to the Lyman continuum photons before the end 
of reionization, every single hydrogen-ionizing photon will be absorbed and converted to 
line emission; thus, the line contribution should be proportional to a production rate of 
hydrogen- ionizing photons, Qh, as 

niijn) = fiQ H (m), (18) 

where /j is a fraction of ionizing photons which are converted to a line i. 

Which lines are important to the near infrared background? The Lyman series photons 
are in right bands; however, they are strongly absorbed and eventually converted to other 
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lines. One exception is the Lyman-o photons: while they are also strongly absorbed, they 
are re-emitted back again in Lyman-o. Therefore the net effect is that the Lyman-o photons 
are not destroyed (in the absence of dust) but merely scattered. Loeb & Rybicki (1999) have 
shown that as the universe expands the Lyman-o; photons are eventually "redshifted out" of 
scattering and escape freely. The Balmer series photons (and others) have too low an energy 
to be relevant to the near infrared background (a direct recombination to n = 2 results in 
a line at 3.4 eV or 3648 A, which will be redshifted to A > 3 /im and is thus irrelevant). 
Therefore, we consider only Lyman-o photons: 



where hv\ yoi = 10.2 eV, fi ya = 0.64, and v Xyoi = 2465 THz. Note that fi ya = 0.64 was derived 
as follows: every hydrogen-ionizing photon results in a n = 2 — > 1 transition. (This is because 
every electron that goes directly to the ground state from n > 3 emits Lyman-series photons 
which are strongly absorbed, creating another excited atom, and this process repeats until 
all electrons end up in n = 2 state.) About 2/3 of the time a n = 2 — > 1 transition creates 
a Lyman-o photon via 2p — > Is transition and about 1/3 of the time it emits continuum 
emission via 2-photon decay of 2s — > Is. The precise value of fi ya depends slightly on the 
temperature of gas, and for a gas at 20,000 K the value of fi ya is 0.64 (Spitzer 1978). Finally, 
we ignore helium recombination lines as their flux is at most 6% of the hydrogen-ionizing 
flux even for metal- free stars (see Table 1 and 4 of Schaerer (2002)). As for a line profile, we 
take it to be a delta function: 



This is an excellent approximation as we are interested in the background intensity which 
is integrated over a broad range of redshifts. If we are, on the other hand, interested in 
a spectrum of individual objects with fine spectral resolution, more accurate calculations 
are required (Loeb & Rybicki 1999). We have confirmed validity of our approximation by 
comparing the resulting spectrum with and without the exact line profile taken into account. 

It should be emphasized that the escape fraction, a fraction of ionizing photons escaping 
from nebula, does not alter luminosity of Lyman-o: very much. This is because all of the 
ionizing photons will eventually be converted to Lyman-o photons which, in turn, will escape 
freely via the cosmological redshift. Therefore, our prediction is free from uncertainty in the 
escape fraction. In other words, we do not care where those Lyman-o photons come from as 
far as the mean intensity is concerned 3 . If most of the ionizing photons escape from nebulae, 




= fi ya hv lya 4)(v - vi ya )Q H (m) 
~ 4.51 x 10 22 erg s" 1 Hz" 1 - 




(19) 




(20) 



3 However, the escape fraction should affect fluctuations as it changes morphology of the ionized region. 
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Lyman-a photons should come from the intergalactic medium (but not too far away from 
nebulae; otherwise, the Lyman-a signal would be spread over a large frequency range and 
the signal would be suppressed). If none of the ionizing photons escape, Lyman-a photons 
should come from nebulae. In both cases, the resulting flux in Lyman-a should be about 
the same. 



2.5. Two-photon Emission 

Luminosity of 2-photon emission is given by 



L v (m) 



Ihv ^ 



fi ya )P{v/vi ya )Q H {m) 



~ 4.24 x 10 22 erg s" 1 Hz" 1 



Qh{ 



in 



10 49 s- 



Vlya/2 



(21) 



where P(y)dy is the normalized probability per 2-photon decay of getting one photon in the 
range dy = du/ui ya . This formula is easily understood: again, every single ionizing photon 
results in a n = 2 — > 1 transition, and I — fi ya — 1/3 (more precisely 0.36 for T g = 20, 000 K) 
of time it emits 2 photons via 2-photon decay. (Therefore there is a factor of 2 multiplying 
his.) We have fitted the data given in Table 4 of Brown & Mathews (1970) to obtain 4 



P(y) = 1.307 - 2.627(y - 0.5) 2 + 2.563(y - 0.5) 4 - 51.69(y - 0.5) 6 . 



(22) 



Note that P(y) is normalized such that Jjj 1 P(y)dy = 1. (This fitting formula gives P(y)dy 
1.0047.) 



3. ENERGY BUDGET 

3.1. Initial Mass Spectrum 

In order to calculate a typical spectrum of radiative efficiency, (e„) (Eq. [3]), one needs 
to specify the mass spectrum of stars, f(m), which determines the mean stellar mass of star 
formation. (For the precise definition of f(m), see Appendix A.) This is important because, 
depending on which mass is the most typical one, hardness of the emerging stellar spectrum 
changes significantly. (Hardness affects the ratio of energy in Lyman-a to that in stellar 

4 Brown & Mathews (1970) tabulate g(y) = 2hyP(y), where h = 6.626 x 1CT 27 erg Hz" 1 is Planck's 
constant. 
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continuum.) Unfortunately, since we have very little knowledge of f{m) for early stars, we 
are not able to estimate the proper error that would result in changing f(m). While we try 
to explore a range of models for f(m), one should keep in mind that our exploration of the 
form of /(to) is not exhaustive. 

We use three different mass spectra: (a) Salpeter (Salpeter 1955): 

f(m) oc m~ 2 - 35 , (23) 

and (b) Larson (Larson 1998): 

/ \ -1.35 

f(m) oc m- 1 1 + — , (24) 

\ m c / 

which matches Salpeter's in the limit of m c — > 0, and one can explore a variety of models 

by changing one parameter, m c . We shall assume m c = 50 M . Finally, (c) a top-heavy 
spectrum: 

> f m-\ 100 < m < 500 M , n . 

t(m) oc < , 25 

J y ' \ 0, otherwise v ; 

which might be possible for the primordial metal-free stars (Bromm & Larson 2004). (Note 
that mf(m) is flat for 100 < m < 500 M Q .) The normalizations are given by 

/ dm f{m) = 1. (26) 

J mi 

The choice of the mass range is somewhat arbitrary. Throughout this paper, we shall assume 
m\ = 5 M and m 2 = 150 M for the Salpeter and Larson mass spectra, whereas m± = 
100 M and m2 = 500 M for the top-heavy spectrum. (We shall explain the reason for 
mi = 5 M later in § 4.2.) The mean stellar masses (Eq. [A5]) are 13.6, 27.4, and 248.5 M 
for the Salpeter, Larson, and top-heavy spectrum, respectively. 



3.2. Metallicity 

The next ingredients are the stellar luminosity- mass relation, L*(m), the ionizing flux- 
mass relation, Q H (m), the stellar lifetime- mass relation, r(m), and the effective temperature- 
mass relation, T e ff(m). Since these relations mainly depend on metallicity, we explore two 
cases: (1) metal-free (Z = 0) stars, and (2) metal-poor (Z = 1/50Z Q ) stars. 
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Fig. 1. — Stellar properties. (Top-left) The stellar bolometric luminosity, Ll ol ; (Top-right) 
the stellar main-sequence lifetime, r; (Middle-left) the number of hydrogen-ionizing photons 
emitted per unit time, Q H ; (Middle-right) the stellar effective temperature, T e ff] (Bottom- 
left) nuclear burning efficiency, with the axis on the right showing bolometric radiative energy 
per stellar mass over the main sequence lifetime; and (Bottom-right) the total number of 
hydrogen-ionizing photons per unit stellar mass, for different stellar metallicities. In the 
top-right panel, the dashed and dotted horizontal lines show the age of the universe at z = 7 
and 15, respectively. 
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3.2.1. Metal- free Stars 

Schaerer (2002) has calculated emission properties of metal-free stars. We use his fitting 
formulas for Q H (m) and r(m) over m = 5 — 500 M & (Table 6 of Schaerer (2002)): 

, n f 43.61 + 4.90a; - 0.83a; 2 9 - 500 M , 

log in \Qh h 1 = S (27) 
&1 ° LV ^ 7 J \ 39.29 + 8.55a; 5 - 9 M , v ; 

log 10 [r/yr] = 9.785 - 3.759a; + 1.413a; 2 - 0.186a; 3 , (28) 
where x = log 10 (m/M e ). We calculate the stellar radius as 

where a = 5.67 x 10~ 5 erg s^ 1 cm~ 2 K+ 4 is the Stephan-Boltzmann constant. The bolometric 
stellar luminosity before absorption, L\ ol , and the effective temperature, T e ff, are given in 
Table 3 of Schaerer (2002). These were used to obtain fitting formulas for T e ff aud L^ ol , 
which are good for masses anywhere from 5 — 1000 M : 

log 10 [L* bol /L e ] = 0.4568 + 3.897a; - 0.5297a; 2 (30) 
log 10 [T e ff/K] = 3.639 + 1.501a; - 0.5561a; 2 + 0.07005a; 3 (31) 

Note that L\ ol and T e ff were calculated for the zero- age main sequence stars, whereas Q H 
has been averaged over the main-sequence lifetime. Strictly speaking, the former quantities 
should have also been averaged over the stellar lifetime; however, we shall ignore such a 
subtlety and use the zero-age values as tabulated in Schaerer (2002). 



3.2.2. Metal-poor Stars 

For stars with Z = 1/50Z Q , we use the fitting formula for Q H (m) and r(m) given in 
Table 6 of Schaerer (2002): 

logic = 27-80 + 30.68a; - 14.80a; 2 + 2.50a; 3 , (32) 

log 10 [r/yr] = 9.59 - 2.79a; + 0.63a; 2 , (33) 

where x = log 10 (m/M Q ). The formula for stellar lifetimes and ionizing photons is good from 
7— 150 M ; we shall extrapolate it down to 5 M . We calculate the stellar radius (Eq. [29]) 
using Ll ol {m) and T e ff(m). This was fit from stellar models given in Lejeune & Schaerer 
(2001). The fitting formulas were thus obtained as 

log lo [L^/L ] = 0.138 + 4.28a; - 0.653a; 2 , (34) 
log 10 [T e// /K] = 3.92 + 0.704a; -0.138a; 2 . (35) 
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Again, L* bol and T e ff were calculated for the zero-age main sequence stars, whereas Q H has 
been averaged over the main-sequence lifetime. 

3.2.3. Stellar Properties 

Figure 1 shows the bolometric stellar luminosity before absorption, L\ ol (top- left panel), 
the main-sequence lifetime, r (top-right), the number of hydrogen-ionizing photons per sec- 
ond, Q H (middle-left), and the stellar effective temperature, T e ff (middle-right), for Z = 
(labeled as "metal-free") and Z = 1/50Z & ("metal-poor"). The bolometric luminosity is very 
similar for metal-free and metal-poor stars at the same stellar mass down to ~ 10 M Q (Tum- 
linson & Shull 2000), and is almost identical for more massive stars (> 100 M Q ) (Bromm, 
Kudritzki & Loeb 2001; Abel, Bryan & Norman 2002). Since metal-free stars had to begin 
their nuclear burning via the p-p chain, which is less efficient than the CNO cycle because 
of weak interactions, the temperature of metal-free stars must be maintained higher than 
that of metal-poor stars to prevent gravitational collapse (Tumlinson & Shull 2000). Since 
the luminosity is similar, this property makes the size of metal-free stars smaller and the 
main-sequence lifetime slightly shorter than those of metal-poor stars. On the other hand, 
metal-free stars produce more hydrogen-ionizing photons than metal-poor stars, particularly 
for m < 30 M Q , owing to their higher temperature (the spectrum is harder). 

The bottom panels of Figure 1 show quantities more relevant to the radiative efficiency, 
e u (m). The first panel shows the ratio of the stellar bolometric energy to the rest mass energy. 
This figure shows that for m ~ 10 — 100 M anywhere from 0.1 to 0.3% of the rest mass 
energy of the star goes into radiative energy via nuclear fusion; thus, this quantity represent 
a "nuclear burning efficiency" of stars. The metal-poor stars radiate slightly more energy 
over their lifetime than the metal-free stars, as they live slightly longer and the bolometric 
luminosity is about the same. On the right, the total number of ionizing photons per unit 
stellar mass, Q H r/m, is shown. The metal-poor stars emit significantly less ionizing photons 
for m < 30 M Q : this property becomes important when we interpret the predicted spectrum 
of the near infrared background. 

3.3. Energy Spectrum 

Using these fitting formulas and the initial mass spectra, we calculate a spectrum of 
radiative efficiency averaged over the mass spectrum for m = mi — m 2 (Eq. [3]). Figure 2 
shows v{e") for the stellar (Eq. [6]), nebular continuum (free-free and free-bound) (Eq. [12]), 
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Fig. 2. — Spectrum of radiative efficiency, ^(e^), for different metallicities: Z = (left) 
and Z = 1/50Z Q (right). The radiative efficiency is defined as a fraction of the stellar 
rest mass energy that goes into free- free (the lowest purple lines), free-bound (the "saw- 
shaped" cyan lines), stellar black-body (green), two-photon emission (orange) or Lyman-o; 
line (the vertical red lines) emission. The solid, dot-dashed, and triple-dot-dashed lines show 
the Salpeter (Eq. [23]), Larson (Eq. [24]), and top-heavy (Eq. [25]) initial mass spectra, 
respectively. 
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Lyman-o; (Eq. [19]), and two-photon (Eq. [21]) emission. The nebular continuum dominates 
at low energy, while the stellar, Lyman-a and two-photon emission dominate at high energy, 
as expected from their spectral shape. Since metallicity changes hardness of the stellar 
spectrum, it affects the ratio of energy in Lyman-a and two-photon emission to stellar 
emission energy: the harder the spectrum is, the more the ionizing photons are emitted, 
and thus the more the Lyman-a and two-photon emission emerge 5 . This explains why the 
metal-free stars have much more energy in Lyman-a and two-photon emission than in stellar 
emission. On the other hand, the metal-poor stars have more energy in stellar emission. 
For the same reason, heavier mass spectra tend to produce more energy in Lyman-a and 
two-photon emission than in stellar emission. In both cases, however, the total radiative 
efficiency is about the same: v(e v ) ~ 10~ 3 . This is merely an approximate conservation of 
energy: initially all the energy was generated by nuclear burning in stars. The generated 
energy is then radiated or reprocessed, but the sum should be more or less the same as 
the input energy. (Of course conservation cannot be exact because we have ignored other 
emission processes such as Balmer lines, helium or metal lines, etc. If the HII region expands, 
additional energy would be lost to expansion.) This property makes the prediction of the 
near infrared background very robust, up to an unknown star formation rate, p*, which will 
be constrained by a comparison to the observational data. 



4. SPECTRUM OF THE NEAR INFRARED BACKGROUND 

4.1. Dependence on Metallicity and Initial Mass Spectrum 

By integrating the volume emissivity over redshift, we obtain the background intensity 
spectrum of the near infrared from early stars (Eq. [4]). To do this, however, one needs to 
specify the evolution of star formation rate over time, p*{z), which is unknown. Therefore, 
for simplicity, we shall assume that the star formation rate is constant over time, at least 
for the redshift range of interest. In other words, we calculate the intensity spectrum for a 
given "time-averaged" star formation rate. Figure 3 and 4 show vl v j p* for stars in three 
redshift ranges, z = 7 — 15, 15 — 30, and 7 — 30. These figures clearly show that the intensity 
at 1 — 2 pm is almost entirely determined by the contribution at z = 7 — 15. (Note that 
Lyman-a lines at z = 7 — 15 are redshifted to 1 — 2 /im.) Therefore, the spectrum of the 



5 Howcvcr, this is not always the case. The bottom-right panel of Figure 1 shows that metal-poor stars 
actually emit as many ionizing photons per stellar mass as metal- free stars for m > 30 Mq ; thus, if the mean 
stellar mass of metal-poor stars were 3> 30 M Q , metal-poor stars would result in as many Lyman-a photons 
as metal-free stars. 
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near infrared background at 1 — 2 /im constrains the star formation rate at z = 7 — 15! 

Table 1 summarizes values of vl v j p\ averaged over 1 — 2 /im. Within 1 — 2 /im, the 
intensity is dominated by Lyman-a emission. For metal poor stars, there is also a significant 
contribution from stars themselves, which brings the overall intensity for metal poor and 
metal free stars to be about the same. This seems striking, but is merely a consequence of 
an approximate energy conservation, as discussed in § 3 6 . Therefore, the predicted intensity 
is not sensitive to stellar metallicity. 

As for dependence on the initial mass spectrum, f(m), heavier mass spectra tend to give 
higher background intensities. Energetics implies that dependence of vl v jp\ on metallicity 
or f(m) should be essentially given by that of the nuclear burning efficiency averaged over 
a population of stars. The last column of Table 1 shows the mass-weighted mean nuclear 
burning efficiency, (L^r/ (mc 2 )), which is tightly correlated with the total signal. Therefore, 
one can explore dependence of the near infrared background on these parameters by simply 
calculating the nuclear burning efficiency dependence on these parameters. In order to 
illustrate how nuclear burning efficiency changes with respect to the shape of f(m), we 
show the efficiency for various values of the lower mass limit, m 1; and the critical mass, 
m c , for the Salpeter (Eq. [23]) and Larson (Eq. [24]) initial mass spectra in Figure 5. The 
average nuclear burning efficiency for mi > 20 M & depends very weakly on m c , while the 
dependence is stronger for mi < 20 M & . Dependence on mi also becomes stronger as mi 
becomes smaller. Overall, for Larson's mass spectrum, different mi and m c may change the 
predicted intensity by a factor of a few, but not much more. However, one should keep in 
mind that other forms of f(m) that we have not explored here might change the predicted 
intensity by a larger factor. 

Mass loss from stars can also affect the total spectrum. Schaerer (2002) gives properties 
of metal free stars undergoing strong mass loss. Using these values, nuclear burning efficiency 
would be increased by a factor of 1.25. The strong mass loss, if any, therefore reduces the 
inferred star formation rate by roughly 25-50 %. (§ 5.1) 

4.2. Comparison with Previous Work 

Santos, Bromm, & Kamionkowski (2002) calculated the near infrared background from 
metal-free stars, assuming that all early stars contributing to the background light have 



6 We thank Paul R. Shapiro for pointing out potential importance of metal-poor stars for the near infrared 
background. 
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Fig. 3. — Spectrum of the near infrared background, vl v /p\, for star formation from z = 
7—15 with various assumptions about metallicity and initial mass spectrum. ("SFR" in 
the axis label denotes /?*.) The left panel shows the metal-free case (Z = 0), while the right 
panel shows the metal-poor case (Z = 1/50Z Q ). The solid, dashed, and triple-dotted-dashed 
lines represent the Salpeter, Larson, and top-heavy mass spectrum, respectively. The thick 
black lines show the total spectrum, while the thin purple, cyan, green, red and orange lines 
show individual contribution from free-free, free-bound, stellar, Lyman-a and two-photon 
emission, respectively. 
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Fig. 4. — The same as Figure 3 for z = 15 — 30 (top) and z = 7 — 30 (bottom). 
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Metallicity 


IMF 


Stellar 


Free- free 


Free-bound 


Ly-a 


2-photon 


Total 


10 3 (^H r ) 


Metal-Poor 


Salpeter 


1.46 


0.00761 


0.0590 


2.10 


0.530 


4.16 


1.35 




Larson 


1.60 


0.0141 


0.109 


3.90 


0.984 


6.61 


1.84 


Metal-Free 


Salpeter 


0.678 


0.00979 


0.0759 


2.71 


0.683 


4.15 


1.25 




Larson 


0.605 


0.0154 


0.120 


4.27 


1.08 


6.08 


1.73 




Top Heavy 


0.663 


0.0205 


0.159 


5.67 


1.43 


7.95 


3.23 



Table 1: Values of vl v j p* averaged over 1 — 2 pm in units of nW m -2 sr~ 1 /M Q yr _1 Mpc~ 3 
from each radiative process. Here p* is the mean star formation rate for m > 5 M & at 
z = 7 — 15, and {L* hol r / (mc 2 )) is the mass-weighted mean nuclear burning efficiency of stars. 



E 




20 40 60 80 100 20 40 60 80 100 

m , m 1 



Fig. 5. — Mass-weighted mean nuclear burning efficiency, (L\ o1 t j (mc 2 )), resulting from vary- 
ing the initial mass spectrum. Each contour shows a different nuclear burning efficiency, and 
m x and m c are the lower mass limit of star formation and the critical mass in Larson's mass 
spectrum (Eq. [24]), respectively. 

m = 1000 M Q (i.e., f(m) = 5 D (m - 1000 M )). While the stellar, Lyman-a, and free- 
free emission were included in their calculations, the two-photon and free-bound emission 
were ignored. As we have shown in Figure 3, the contribution from free-bound transition to 
n = 2 is much larger than the free-free contribution, and the contribution from two-photon 
emission is as important as that from stellar emission; thus, they must be included. 
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Salvaterra & Ferrara (2003) assumed metal-free stars but explored different initial mass 
spectra, varying m c of Larson's spectrum. There is a subtle but important difference between 
their approach and our approach. When they modeled the volume emissivity, they did not 
allow for the faster rate of death of higher mass stars. In other words, their formula for 
the volume emissivity implicitly assumed that stars lived longer than the Hubble time. We 
find that the volume emissivity that is not corrected for dead stars (Eq. [A6]) agrees with 
their formula (Eq. [5] and [12] in Salvaterra & Ferrara (2003)). This assumption leads to 
an overestimation of contributions from higher mass stars which live much shorter than the 
Hubble time. The physical reason is because higher mass stars live shorter and die sooner, 
and therefore there are fewer massive stars around to emit energy at a given z. By properly 
taking into account the faster rate of death of stars, we have obtained the correct formula 
(Eq. [A9]) which explicitly depends on the stellar main-sequence lifetime. On the other hand, 
if stars live longer than or comparable to the Hubble time, then our approximation breaks 
down and one should do the integral in Eq. [A8]. This is precisely why we have restricted 
our attention only to fairly massive stars, m > 5 M , for which the age of the universe is 
definitely longer than the main-sequence lifetime at redshifts of interest. 



5. IMPLICATIONS FOR THE COSMIC STAR FORMATION RATE 

5.1. Star Formation Rate at z = 7 — 15 

Comparing the predicted values of vl v jp i( (Table 1) to the measured data, we can 
constrain the star formation rate p*. The near infrared background has been determined 
with various satellites, such as the Diffuse Infrared Background Experiment (DIRBE) on 
the Cosmic Background Explorer (Hauser et al. 1998; Boggess et al. 1992) and the Near 
Infrared Spectrometer (NIRS) on the InfraRed Telescope in Space (IRTS) (Matsumoto et 
al. 2005). Table 2 summarizes the observational data. A significant uncertainty exists in the 
observational data, largely because of uncertainty in subtraction of the zodiacal emission. A 
large difference between Wright (2001) and Cambresy et al. (2001), which have used the same 
data (DIRBE), is entirely due to difference in the zodiacal light models. One may summarize 
the current measurement of the cosmic near infrared background as ~ 2 — 50 nW m -2 sr~ x 
in 1-2 pm, which includes the 1-a lower bound of the lowest measurement and the 1-er 
upper bound of the highest measurement. 7 . Taking into account a scatter in theoretical 



7 A recent analysis of blazar spectrum by Aharonian et al . (2005) suggested that the intensity of the near 
infrared background may be lower, with an upper limit of ~ (14 ± 4) nW m~ 2 sr _1 from 1-2 /im, which 
is more consistent with the analysis of the DIRBE data by Wright (2001). On the other hand, a recent 
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Instrument 


vl v 


Reference 


DIRBE 
NIRS 


14.6 ± 12.4 
30.9 ± 12.8 
41.6 ±9.7 


Wright (2001) 
Cambresy et al. (2001) 
Matsumoto et al. (2005) 



Table 2: Observational data of vl v averaged over 1 — 2 /im in units of nW m~ 2 sr" 1 . Note 
that zodiacal emission, emission of foreground stars, and emission from galaxies have been 
subtracted. For DIRBE, the straight average of J- and K-band data (minus 10 nW mT 2 sr -1 
which is contributed by galaxies) is quoted, and the error is estimated as a/ (<j 2 + a 2 K )/2. 
For NIRS, the straight average of 9 measurements from 1.43 to 2.24 /im is quoted, and the 

error is estimated similarly as \JY1^=i °f/9- Statistical errors are negligible compared to the 
systematic errors in subtraction of the zodiacal emission. For details, see Matsumoto et al. 
(2005). 

predictions due to different assumptions about metallicity and initial mass spectrum (see 
Table 1 and 3), we obtain p* ~ 0.3 - 12 M yr" 1 Mpc" 3 at z = 7 - 15. (Note that the error 
bar is not dominated by theory but by observational errors.) What does this imply? 



5.2. Stellar Mass Density Confronts Cosmic Baryon Density 



One must always make sure that the stellar mass density inferred from star formation 
rate does not exceed the cosmic mean baryon density. Using the formula derived in Ap- 
pendix B, we obtain the ratio of cumulative mass density (which is not corrected for dead 
stars) of stars formed at 7 < z < 15 to the cosmic mean baryon density as 



pf m (7 < z < 15) 



~ 0.0747 



PbO 



0.024 
Q h h 2 



0.14 
Vt m h 



1/2 r 



p«(7 < z < 15) 
M yr -1 Mpc" 3 



(36) 



where pbo = 2.775 x 10 11 (O^/i 2 ) M Q Mpc -3 is the present-day mean baryon density. (Note 
that p* denotes comoving mass density.) It follows from this equation that the inferred lower 
limit to the star formation rate, p* > 0.3 M Q yr" 1 Mpc -3 , requires that more than 2% of 
baryons in the universe should have been converted into stars. Madau & Silk (2005) argue 
that "this is energetically and astrophysically daunting" . It would be daunting, if the stars 
that were responsible for producing the near infrared background lived longer than the age 
of the universe, and more than 2% of baryons had remained locked up in the stars. However, 



detection of the fluctuations in the near infrared background may imply that at least some of the NIRB is 
from early stars (Kashlinsky et al. 2005) 
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Metallicity 


IMF 


p\ (Lower 


p\ (Upper 


Ngen 


P*/Pb [%] (Lower 


P*/Pb [%] (Upper 






I — a limit) 


l — o" limit) 




1 — o limit) 


1 — o limit) 


Metal-Poor 


Salpeter 


0.48 


12 


7.3 


0.49 


12.3 




Larson 


0.30 


7.6 


12.8 


0.18 


4.41 


Metal-Free 


Salpeter 


0.48 


12 


10.0 


0.36 


9.0 




Larson 


0.33 


8.2 


16.7 


0.15 


3.7 




Top Heavy 


0.25 


6.3 


120 


0.016 


0.39 



Table 3: Implications of the current lower and upper 1 — o limits to the cosmic near infrared 
background data, vl v = 2 nW m -2 sr" 1 and vl v = 50 nW m -2 sr _1 , for formation of early 
stars with m > 5 M Q during z = 7 — 15. The third and forth columns show the lower and 
upper limits, respectively, to the star formation rate in units of M yr _1 Mpc~ 3 , the fifth 
column shows the mean number of generations of stars (Eq. [37]), and the sixth and seventh 
columns shows the lower and upper limits, respectively, to a ratio of stellar density to the 
mean cosmic baryon density in percentage. 



it is certain that stars lived much shorter than the age of the universe (see the top-right 
panel of Figure 1), and the stellar mass density must be corrected for dead stars; thus, the 
actual amount of baryons locked up in stars at any given time between z = 7 and 15 should 
be less than that is given by Eq. [36]. We derive the stellar mass density corrected for dead 
stars in Appendix B.2, which shows that the correct answer should lie between Eq. [36] and 
Eq. [36] divided by the mean number of generations of stars, N gen , given by 

N = *( 7 <*< 15 ) (37) 

and t(7 < z < 15) = 266 million years is the cosmic time elapsed during z = 7 — 15. Table 3 
tabulates N gen for various assumptions about metallicity and initial mass spectrum. From 
this we conclude that, to explain the cosmic near infrared background by early generations 
of stars, 0.016-12% of baryons need to be processed in stars at a given time between z = 7 
and 15. If we take the lower l-o limit, only 0.016-0.49% of baryons need to be processed in 
stars (depending on metallicity and mass spectrum); this is not a daunting requirement and 
does not exclude the stellar origin of the cosmic near infrared background 8 . 



8 Our argument so far has implicitly assumed that all of baryonic gas in the previous generation of stars is 
returned to the intergalactic medium and recycled in the subsequent generation of stars. In reality, however, 
only a fraction of gas would be returned (and the rest of gas would be locked up in compact remnants such 
as black holes); thus, the real requirement would be somewhat larger than 0.016-0.49%. 
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5.3. Comparison with Low Redshift Data 

How does the inferred star formation rate at z = 7 — 15 compare to the \ow-z rate? 
Figure 6 compares the cosmic star formation rate at z < 5 (Gabasch et al. 2004) 9 to that 
constrained by the near infrared background. While uncertainty due to subtraction of the 
zodiacal light is large, it is quite clear that the star formation rate at z = 7 — 15 required to 
account for the cosmic near infrared background data is much higher than that at z < 5 by 
more than an order of magnitude. 

It must be emphasized, however, that Figure 6 is potentially misleading: as we have 
already discussed, the star formation rate inferred from the near infrared background is 
only for stars more massive than 5 M . On the other hand, the low-z data are primarily 
dominated by low mass stars; thus, Figure 6 might be comparing apples and oranges. As 
such low mass stars do not contribute to the near infrared background, it is not possible to 
infer their formation rate directly. One may still estimate it by extrapolating the initial mass 
spectrum down to lower masses, and by doing so the total star formation rate at z = 7 — 15 
should rise. In other words, the constraint shown in Figure 6 should be taken as a lower 
bound. Also, dust extinction (which we have ignored), if any, would make the required star 
formation rate rise even higher. 

6. METALLICITY CONSTRAINTS ON STAR FORMATION 

One of the ways to constrain early star formation is to take into account the amount 
of metals that can be produced without over-polluting the universe. Metals ejected from a 
star have two origins: 1. stellar winds, which inject metals into the IGM over the course of 
the star's lifetime; and 2. the final disruption of the star. Stars of low metallicity end their 
lives in different ways and produce different amounts of metals, according to the initial mass 
of the star (Heger et al. 2003; Portinari, Chiosi, & Bressan 1998; Siess, Livio, & Lattanzio 
2002). 

The metal yields of stars with initial metallicity of Z = 1/50Z Q were given in Portinari, 
Chiosi, & Bressan (1998). These models of metal production take into account stellar winds 
and supernova explosions. The metal production efficiency (metal mass ejected from the 
star over initial stellar mass) is shown in Figure 7. It is clear that metal production depends 
strongly on the initial mass of the star and how the star ended its life. 



9 The rate at z < 5 has been shifted upward by 0.35 dex to correct for dust extinction. More recent 
determination of the star formation rate by Drory et al. (2005) agrees very well with Gabasch et al. (2004). 
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Fig. 6. — Cosmic star formation rate. The shaded region shows the star formation rate 
for m > 5 M Q constrained by the current data of the cosmic near infrared background, 
vl v ~ 2 — 50 nW m~ 2 sr -1 , in 1-2 /im. The solid line shows the star formation rate at z < 5 
(Gabasch et al. 2004). Note that the shaded region should be taken as a lower bound. (See 
discussion in § 5.3). 



• From 6 to 8 M & , the O/Ne/Mg core of the star collapses, or the star ejects its outer 
envelope, leaving a white dwarf or neutron star. 

• From 8 to 25 M Q , the iron core collapses, the star explodes as a supernova, and a 
neutron star is left as a remnant. A significant amount of metals are ejected. 

• From 25 to 40 M Q , there is a weak supernova and a black hole is created by fallback. 
The amount of metals that are ejected into the IGM decreases sharply, leaving most 
of the metals locked in the black hole. 

• From 40 to 100 M Q , the star directly collapses into a black hole. The only metals 
produced are from mass loss during the star's life. 

• From 100 to 140 M , a pulsational pair instability supernova results. This ejects the 
outer envelope of the star, and then the core collapses into a black hole. Metals in the 
outer envelope pollute the IGM. 
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• From 140 to 260 M , a pair instability supernova results, which completely disrupts 
the star and leaves no remnant. All the metals are ejected into the IGM. 

• Above 260 M , the star collapses directly into a black hole, and there is no enrichment 
of the IGM. 



Fitting functions were made for each of these mass intervals. 



M meta i S;e j(m) 



m/M G 







0, < 8 M 

-0.0708 + 0.0142^, 8-25M 

0.142, 25 - 40 M 

0.0255, 40 - 100 M ^ 

0.0117 + 0.00152-^, 100 - 140 M 

0.485 - 0.000144^, 140 - 260 M 

M meta is,ej is the mass in solar masses of ejected metals. These fitting functions are limited 
because of the small number of data points being fit. However, they should be able to serve 
to approximate the amount of metals formed in each mass range. For stars less massive than 
8M , the amount of metals is negligible, and is assumed to be zero. When pairing these 
metal production rates with an IMF, they agree well with Giroux & Shapiro (1996). 

Combining this with an initial mass function, we are in the position to be able to predict 
how many metals are formed by these first stars in a certain redshift range. The amount of 
metals formed from stars is given by 

Pmetais = — / / (m) M metals>ej (m)dm (39) 

where p me tais is the mass density of metals from a population of stars, p* is the stellar mass 
density, and is the mean stellar mass of the mass function (Eq. [A5]). 

The very first stars were metal free. We will assume that metal free star formation 
only occurs to a certain critical metallicity, Z = 10~ 4 Z & . (Schneider et al. 2002; Bromm 
& Larson 2004). After this metallicity is reached, metal poor stars begin to form. If this 
transition happened before a redshift of 15, all of the NIRB is from metal poor stars, and the 
maximum amount of metals would be produced. Would these metal poor stars overproduce 
metals that we observe in the universe today? 

We take the upper limit of metals that can be produced in the universe to be the solar 
metallicity, Z = Z & = 0.02. Using Q b = 0.044 and H = 70 kms _1 Mpc _1 , the baryonic 
density, pb, is 6.0 x 10 9 M Mpc~ 3 . Then the maximum metal mass density that can be 
produced is then 

Petals = P b Z Q (40) 
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Fig. 7. — The ratio of ejected stellar metals and stellar mass versus the mass of the star. 
Note that each jump corresponds to a different stellar fate. The lines correspond to the 
fitting functions given in Eq [38]. Open diamonds correspond to the stellar models given in 
Portinari, Chiosi, & Bressan (1998). 
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IMF 


P* 


Pmetals M Q MpC 3 


Salpeter 


0.48 


9.59 x 10 5 




12 


2.40 x 10 7 


Larson 


0.30 


1.61 x 10 7 




7.6 


4.09 x 10 s 



Table 4: The mass density of metals being produced by metal poor stars from a redshift of 
7 to 15. Values of p\ are the upper and lower 1-a limits given in Table 3. All cases (except 
for the Larson mass function with an upper 1-a limit for p\) give values of Pmetals that is less 
that the metal density we see today (p me tais — 1-2 x 10 8 M Q Mpc~ 3 ), indicating that metal 
poor stars can form up to z = 7 and not overproduce metals. 

This gives a metal mass density of petals = 1-2 x 10 s M Mpc~ 3 . If the predicted metallicity 
from Eq. [39] exceeds this value, the model must be ruled out. 

A population of metal poor stars that form from z=7-15 will form metals according to 
Eq. [39]. This depends on the mass function and the star formation rate. The amount of 
metals formed for the upper and lower 1 — a limit of p\ for each mass function are shown 
in table 4. We have found that a population of metal poor stars do not overproduce metals 
that we observe in the universe today, except for the Larson mass function upper 1 — a value 
for the star formation rate. 

In a recent paper by Dwek, Arendt & Krennrich (2005), they showed if the entire NIRB 
excess is from stars, these stars would produce about 87 % of the metal and helium abundance 
we see today. This would leave little room for production in more recent eras from Pop II 
and Pop I stars. In order to avoid this, they stated that stars can collapse directly into black 
holes, locking up their metals and therefore not ejecting any metals into the IGM. They 
used this arguement to suggest that the NIRB must come from very massive Pop III stars, 
which are massive enough to collapse directly into black holes. We showed that this is not 
necessarily the case: metal poor stars with a regular mass function can also produce less 
than the observed amount of metals. By doing the calculations using the stellar data that 
takes into account both ejected metals and remnants, we find that the first stars may only 
contribute as little as less than 1 % of the metals in the IGM, which is certainly feesable. 

We now have a picture of how star formation could have occurred in the early universe. 
Metal free stars began forming early. As they die, they pollute the IGM with metals. Once 
a critical metallicity is reached, metal poor stars begin to form. These metal poor stars can 
continue to form up to a redshift of 7 without overproducing metals in the IGM, while still 
providing enough ultraviolet photons to make up the Near Infrared Background. 
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7. AMPLITUDE OF MATTER FLUCTUATIONS AT SMALL SCALES 

In the previous work on the near infrared background from early stars (Santos, Bromm, 
& Kamionkowski 2002; Salvaterra & Ferrara 2003; Madau & Silk 2005), it was concluded 
that a substantial fraction of collapsed baryons in the universe must have been converted into 
stars. In § 5.2 we have considered a fraction of all baryons in the universe that was locked 
up in stars at a given z, and now it is natural for us to ask how many baryons collapsed in 
dark matter halos were converted into stars. Of course, this fraction, rj, which is sometimes 
called the "star formation efficiency" , must not exceed unity. 

In Appendix C, we derive an analytical model for the cosmic star formation rate 
(Eq. [C6]) which has often been used in the literature: 

/*.(*) ~ 0.536 M Mpc- 3 yr- 1 (X) (i±£ ^ y min (z) e -^/ 2 , (41) 

where 

V W - (42) 

and <t(M) is the present-day r.m.s. amplitude of mass fluctuations given by 

3ji(kR M ) 



ro../i.n \ "1 2 



kR 



M 



(43) 



Here, P(k) is the linear power spectrum of density fluctuations at present, ji(x) = sin(:r)/x 2 — 
cos(x)/x is the spherical Bessel function of order 1, R M is a radius defined by M = 
(47r/3)p m0 i?5^, and M min (z) is the minimum mass of dark matter halos which can host star 
formation. For illustration purposes, let us assume that star formation occurs when cooling 
via atomic hydrogen becomes efficient, M > M min (z) ~ 10 8 M Q [(1 + z)/10]~ 3//2 (Barkana 
& Loeb 2001). This mass scale roughly corresponds to the wavenumber of k — h/Rm ^ 

By comparing this analytical model to the observational data, one can constrain r] 
and/or cr[M min (z)]. The usual approach is to constrain r] by fixing o"[M min (z)] with P(k) 
extrapolated from k < 1 Mpc -1 (which corresponds to M > 5 x 10 12 M ) down to much 
smaller scales, k > 36.8[(1 + z)/!^] 1 / 2 Mpc" 1 , for which we do not have any direct obser- 
vational constraints yet. This is potentially a dangerous approach. Eq. [41] implies that 
p*(z) is exponentially sensitive to cr[M min (z)], and a slight increase in cr[M m i n (z)] may sub- 
stantially reduce rj for a given p*(z), and thus one must not ignore the fact that we do not 
know the precise value of <r[M min (z)] without relying on extrapolations of P(k) by a factor 
of more than 40 in k. Therefore, the stellar origin of the near infrared background cannot 
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be ruled out on the basis of 77 until we understand the amplitude of matter fluctuations on 
small scales. One may reverse the argument: it might be possible to explore the small-scale 
fluctuations by using the near infrared background, for a given star formation efficiency 

8. CONCLUSIONS 

We have presented detailed theoretical calculations of the intensity of the cosmic near 
infrared background from early stars. We have shown that the intensity is essentially deter- 
mined by the mass-weighted mean nuclear burning energy of stars (for a given mass spectrum 
of early stars) and the cosmic star formation rate. The prediction is not sensitive to stellar 
metallicity (Table 1), while uncertainty from the initial mass spectrum could be large, as 
we have very little knowledge about the form of mass spectrum for early stars. Our simple 
analytical calculations agree well with recent numerical calculations of the spectrum using 
the CLOUDY code (Dwek, Arendt & Krennrich 2005). The measured intensity at 1 — 2 /j,m 
can be used to infer the cosmic star formation rate at 7 < z < 15, which is difficult to 
constrain by other means. Although the current data are quite uncertain due to subtrac- 
tion of the zodiacal light, the inferred star formation rate, p* ~ 0.3 — 12 M Q yr~ x Mpc~ 3 
at z = 7 — 15 for m > 5 M , is significantly higher than the low- z rates. We have shown 
that this does not exclude the stellar origin of the cosmic near infrared background, as it 
merely requires more than 0.016-0.49% of baryons to be processed in stars at any given time 
between z = 7 and 15 (depending on metallicity and initial mass spectrum; see Table 3). 
Such a high star formation rate at high z may be consistent with recent theoretical propos- 
als (Cen 2003; Mackey, Bromm, & Hernquist 2003). In addition, the derived star formation 
rate does not overproduce metals in the IGM (unless using the upper 1-er value paired with 
a Larson mass function), and may produce as little as less than 1 % of the metals in the 
IGM. More accurate determination of the near infrared background is absolutely necessary 
to yield a meaningful estimate of the star formation rate with any confidence. If the future 
data demand too high of a star formation rate for the stellar origin to be viable, then other 
sources that might contribute to the near infrared background, such as early quasars, should 
be invoked (Cooray k Yoshida 2004; Madau & Silk 2005). 

We thank K. Ahn, M. Alvarez, V. Bromm, N.J. Evans II, A. Ferrara, J. Scalo, P.R. 
Shapiro and N. Yoshida for valuable comments on early versions of the paper and fruitful 
discussions. 
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A. DERIVATION OF VOLUME EMISSIVITY 

In this Appendix we derive the formula for volume emissivity, p(u, z), which is formally 
given by 

/ x f m2 i dnjm, z)—a. . /A1 . 

p(v,z) = }2j dm ^ J L u (m), (Al) 

where a labels various radiative processes (e.g., line), L° is a time-averaged energy spectrum 
of relevant emission (luminosity in frequency interval [v, u+du]), and dn*/dm is the comoving 
number density of stars that are alive at a given z in mass interval [to, m + dm] . Subtlety 
exists as we need to properly take into account dead stars. (Since dead stars do not contribute 
to the emissivity, they must be removed from Eq. [Al].) The goal of this Appendix is to 
derive the formula for p(u, z) which is properly corrected for dead stars. 



A.l. Case with No Dead Stars 

First, as the simplest example (and for illustration purposes) let us derive the formula 
which does not correct for dead stars. (Note that we do not use this formula. The correct 
formula will be given in the next subsection.) We write dn*/dm as 

^*=».WH, (A2) 

where /(to) is a probability distribution function of stellar masses (also known as the mass 
spectrum) normalized to unity for a certain mass range such that 

dm /(to) = l. (A3) 



/ 

J m 



We assume that /(to) is independent of time. For example, /(to) oc m~ 2 ' 35 , the Salpeter 
mass spectrum, is independent of time. In principle, however, /(to) may depend on time 
when there is a characteristic stellar mass scale (e.g., Larson's mass spectrum) that increases 
or decreases with time. One may expect m c to decrease as the metal enrichment proceeds, 
for example. Nevertheless, we shall assume that m c is independent of time, at least for the 
redshift range that we consider. 

We may use the comoving mass density of stars, p*(z), instead of the comoving number 
density, n*(z). The relation is 

p*(z) = m m n*(z), (A4) 
where to* is the mean stellar mass given by 

pm.2 

m»E I dm to/(to). (A5) 

J mi 
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Using this, the volume emissivity (not corrected for dead stars) becomes 

pfa z ) = PM J2 r dmf(m)T v {m). (A6) 

m * a J mi 

Salvaterra & Ferrara (2003) used a version of this Eq. 10 , and therefore they did not correct 
their emissivity for dead stars. 



A. 2. Emissivity Corrected for Dead Stars 

Now, we correct emissivity for dead stars by simply removing them from n*{z): 



f(m) I dt n„, (A7) 



dn*(m, z) 

dm JT (z)-T{m) 

where n* is a rate of star formation, T (z) is the time between when the universe started 
forming stars and the time corresponding to z, and r(m) is a stellar main-sequence lifetime. 
Eq. [Al] then becomes 

P( v > z ) = 22 / dm f(m)L a v (m) I -. (A8) 

a J mi JT (z)-T(m) m * 

This result may be simplified when the stellar lifetime is much shorter than T Q (z) (which is 
about the same as the age of the universe at z). As we have shown in Figure 1, the main- 
sequence lifetime of stars contributing to the near infrared background is always shorter than 
the age of the universe. We Taylor expand the integral over time to obtain the final formula: 

p(u, z) ~ ^ ^ dm f{m)T v {m)r{m). (A9) 

m * J mi 



B. STELLAR MASS DENSITY 
B.l. Case with No Dead Stars 

Cumulative mass density of stars that were formed at t\ < t < t 2 is given by 

pT m (h, t 2 )= f 2 dt p* w (t 2 - t x )p„ (Bi) 
Jtl 



10 Salvaterra & Ferrara (2003) actually used p(i>, z) = p*(z) J2 a C, 2 dmf(m)L"(m)/m, which is off by a 
factor of m/m*. 



-32 - 



where we have assumed that p\ is approximately constant over t\ < t < t 2 . Using z instead 
of t, one gets 

dz i>. 



cum ( 

P* [ z l 



J Z\ 



H(z)(l + z) 



1 



dz p\ 



qU 2 h J Zi (i + ^) 5/2 

5.509 x 10 s M Q Mpc - ' 



0.14 

fit, 



1/2 



P* 



M yr" 1 Mpc- 3 



10 



1 + 21 



3/2 



10 



l + 2 2 



3/2' 



(B2) 



where we have assumed that z 1 3> (^a/^) 1 / 3 — 1 ^ 0.3 so that dark energy contribution to 
H(z) can be ignored. This is cumulative density, as it includes those stars which had already 
died. The fractional stellar mass density relative to the critical density, is then given by 



P* (21,22) 

"* \Zi,Z2) = 



flt um {z u z 2 )h' 



PcO 

1.985 x 10~ 3 
10 

x 



0.14 



1+2! 



3/2 



1/2 



10 



p* 



M yr" 1 Mpc- 3 

\ 3/2' 



1 +2 2 



(B3) 
(B4) 



where p c0 = 2.775/i 2 x 10 11 M Q Mpc 3 is the present-day critical density. 



B.2. Stellar Mass Density Corrected for Dead Stars 

We remove dead stars (whose lifetime, r(m), is much shorter than the Hubble time) 
from the stellar mass density to obtain 

P*(ti,t 2 ) ~ p* / dm /(m)r(m). (B5) 

J mi 

Comparing this with the cumulative stellar density (Eq. [Bl]), one finds the relation 

P*(ti,t 2 ) * ^J? 1 '^ , (B6) 

^ gen 

where N gen is the average number of generation of stars: 

Ngen = rm 2 , 77 \ / \ " (B7) 

J dm /(m)r(m) 
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C. ANALYTICAL MODEL OF THE COSMIC STAR FORMATION RATE 

A popular assumption usually made for an analytical model of the cosmic star formation 
rate at high z is that p\{z) is related to the mass function of dark matter halos: 

• ( \ - ui vi ^_ , dF[M > M min (z)] 

P*{z) =r] Pm H(z)(l + z) — , (CI) 

where p b0 = 2.775 x 10 11 [Vt^h 2 ) M Q Mpc" 3 is the present-day mean baryon density, and r\ 
represents a "star formation efficiency", a constant fraction of baryonic gas in dark matter 
halos that was converted into stars. It is admittedly too simplistic to assume that i] is 
independent of halo mass or redshift. For example, negative feedback from a star forming 
in a single mini-halo might prevent the formation of multiple stars in the same halo, which 
would imply rj oc M~ l . Therefore, this parameterization of the efficiency serves merely as an 
order-of-magnitude representation of the true efficiency. Finally, F is a "collapse fraction" , 
a fraction of mass in the universe collapsed into halos more massive than M min (z), 

F[M > M mm (z)} = ± dM^j^M, (C2) 

where p m0 = 2.775 x 10 11 (Q m h 2 ) M & Mpc -3 is the present-day mean total mass density, 
and M is the halo mass (not to be confused with the stellar mass, m), and drih/dM is the 
halo mass function. The Press-Schechter mass function gives the collapse fraction in terms 
of the complimentary error function, 

F[M > M miQ (z)] = J- I" dy e"^/ 2 , (C3) 

where 

< \ = 6 < ~ 1-68(1 + *) 

VmhAZ} ~ a[M miQ (z))D(z) ~ a[M min (z)} ' 1 } 

in the redshift range of interest (z > 7). One thus finds 

dF[M > M min (z)] [2d yminc _ yLjz)/2 
dz V vr dz 



2 yminV^) 
7T 1 + Z 



^ d In M m \ n (z) din a 1 



d\n(l + z) dlnM 



M=M min (z) 



e -ldlnW/2 



2 y mi n(z) 1 - n cS [M m - m (z)] 

7T 1 + Z 4 



e -i£in(*)/2. (C5) 



Here, we have assumed that M min (z) is proportional to the virial mass and thus M min (z) oc 
(1 + z)~ 3 / 2 , and defined an effective slope of the power spectrum as [n eff (M) + 3]/6 = 



34 



d \na~ l / d\nM . At the scale of interest, the effective slope lies between —2.5 and —2.8; 
thus, the prefactor, (1 — n e s)/A, is of order unity and the exact value is not so important. 
Setting it to be unity, we obtain a fully analytical formula for the star formation rate at high 
redshifts, 
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